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Abstract 

Advances in our numerical and theoretical understanding of gamma-ray burst afterglow processes allow us to construct models 
capable of dealing with complex relativistic jet dynamics and non-thermal emission, that can be compared directly to data from 
instruments such as Swift. Because afterglow blast waves and power law spectra are intrinsically scale-invariant under changes 
of explosion energy and medium density, templates can be generated from large-scale hydrodynamics simulations. This allows 
for iterative template-based model fitting using the physical model parameters (quantifying the properties of the burster, emission 
and observer) directly as fit variables. Here I review how such an approach to afterglow analysis works in practice, paying special 
attention to the underlying model assumptions, possibilities, caveats and limitations of this type of analysis. Because some model 
parameters can be degenerate in certain regions of parameter space, or unconstrained if data in a limited number of a bands is 
available, a Bayesian approach is a natural fit. The main features of the standard afterglow model are reviewed in detail. 


1. Introduction 


Gamma-ray bursts (GRBs) are the most luminous explo¬ 
sions occurring in the universe and a key target for many 
active and upcoming rapid cadence survey programs from 
radio to gamma-rays. We know they most likely involve 
the merging of neutron stars (making them prime gravita¬ 
tional wave counterpart candidates) and massive star collapse, 
and that some fo rm of relativi s tic ou t flow is launched. The 
merger scenario ( Eichler et al. . 1989t Paczvnskl 1991 1 has 
been tied to short GRBs, while l ong GRBs are securely con¬ 
necte d to massive star collapse ( Wooslev . 1993t Paczvnski . 
1998 : MacEadven and Wooslevl 1999h thr ough observations 
of coincident supernovae (IGalamal Il998l) . The separation 
between ‘short’ and ‘lon g’ burst durations lies around 2 s. 
dKouveliotou et al. . 1993). 

We also know this blast wave eventually generates an after¬ 
glow from X-rays to radio, as was predicted originally by way 
of the fireball model which explained GRB prompt emission 
in terms of optically thin synchrotron emission from colliding 
shells within a hydrodynamically launched relat i vistic outflow 
dRees and Meszarc^ 19921: Meszaros and Rees . 1997h . Dur¬ 
ing the afterglow stage the emission is dominated by the for¬ 
ward shock interacting with the external medium, with elec¬ 
trons being accelerated to relativistic velocities at the shock 
front and interacting with the locally generated magnetic held 
in order to produce synchrotron emission. Afterglow emis¬ 
sion is not unique to the synchrotron internal shock model, 
and a decelerating afterglow-stage blast wave can be associ¬ 
ated with a number of mechanisms for jet launchin g (such 
as magnetically dominated jets, e.g. Drenkhahn 2002h and/or 
prompt emission (such a s photospheric emission models, e.g. 
Meszkros and ReMbOOd) . 
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Afterglow observations have proven extremely useful for 
a number of reasons and have been inst rumental in estab¬ 
lishing the extra-galacti c nature of GRB dCosta et all 1 1997 


van 


Paradiis et al. . 1997h . They provide insight in the envi¬ 


ronment of the burster (constraining circumburst medium den¬ 
sity, amount of dust-extinction), the physical properties of the 
progenitor (via explosion energy), the physics of jet launch¬ 
ing (via jet collimation angle), and the fundamental plasma ki¬ 
netic theory of relativistic shocks (via micro-physical parame¬ 
ters describing magnetic held generation and electron acceler¬ 
ation). Most of this is done through interpreting the evolution 
of the characteristics of the synchrotron spectrum in terms of 
flux equations derived from dynamica l blast wave models in 
the se lf-similar ultra-relativistic stage ([BlandfordandMcKee, 
1^7^ the la te self-similar non-relativistic stage ( Tavloii 195C ; 


Sedo 


^ Il959l) . 


or, in a more recent development, from complex 


multi-dimensional, trans-relativistic evolution in-between. 

In this review, I describe the basic aspects of afterglow mod¬ 
els, and the challenges involved in performing data analysis 
based directly on physical models (as opposed to post-hoc in¬ 
terpretation of purely heuristic functions such as power laws, 
that capture the shape of the data of light curves and spectra 
in a simplified manner). In section |2l the dynamics of blast 
waves from ultra-relativistic to non-relativistic are discussed 
and self-contained models are provided. Afterglow emission 
is discussed in section |3] and it is reviewed how flux equations 
for afterglow spectra are derived. Hydrodynamical simulations 
and model-based afterglow fitting are discussed in section |4] 


2. The dynamics of afterglow blast waves 

2.7. Simple approximations for spherical flow 

Under the assumption of spherical flow, the radius R of the 
afterglow blast wave is fixed from conservation of explosion en¬ 
ergy E in the blast wave (assuming adiabatic expansion). The 
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conditions at the front of the shock are set by the shock-jump 
conditions. There exists a reasonable approximation to the ideal 
gas equation of state (EOS) that covers the tra nsition between 


relativ istic and non-relativistic temperatures (IMignone et al 

2nn5h : 

el(pc^) 2 + e/{pc^) 


pKpc ) = 


( 1 ) 


1 + ejipc^) ’ 

where p pressure, c speed of light, p co-moving density, e in¬ 
ternal energy density excluding rest mass. Using this EOS, the 
shock-jump conditions can be simplified to 


P 

e 
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( 2 ) 


valid throughout the evolution of the blast wave (and where 7 
the fluid Lorentz factor at the shock front, pe^t the mass den¬ 
sity in front of the shock and E the Lorentz factor of the shock). 
The upstream density can be allowed to depend on radius, i.e 
Pext = Pref{RIRref) with k - (0,2) Covering respectively 
a homogeneous interstellar-medium (ISM) type and a stellar- 
wind type environment. The concis e des c riptio n offered by 
eq. |2] was pointed out previously by lUhml (1201 ih : Ivan Eerten 
(2013h. and this particular EOS has been used in numerical 


(IZhang and MacFadvenI 120091: Ivan Eerten et al.L 1201 Obi 1201 Ih 


and theoretical ( Nava et al. . 2013h analysis. Although in the 


non-relativistic limit the jump in density is fixed (depending on 
the polytropic index; the jump of 4 from the equations above 
applies to an ideal gas), the shock-jump conditions indicate a 
special feature of relativistic flows, where the density jump can 
become arbitrarily high. This effect is even stronger when ex¬ 
pressed in the lab frame, where an additional factor of y applies. 

Further assuming a homogeneous shell model, the radius of 
the blast wave can be expressed exactly as ( van Eert^ 2013h : 


£/(Mc2) = ^02(4^2 _ 


(3) 


where M the swept-up mass (i.e. a proxy for radius) and p the 
fluid velocity v in units of c. The width of the shell will always 
be A/? = Rl{\2y^). This can be established by taking the down¬ 
stream density according to the shock-jump conditions and, un¬ 
der the assumption of a homogeneous shell, moving down in 
radius until all swept-up matter is accounted for. Due to its in¬ 
clusion of a simplified EOS, the blast wave model described by 
eq. a is about the most concise analytically tractable approxi¬ 
mation possible. 

Other simplified trans- relativistic spherical shell mode ls ex¬ 
ist in the literatu re (e.g. Piran 1999t Huang et al. 19991 who 
omit pressure, or 


Pe’er 


20121 who includes pressure); in prac¬ 


tice, the differences between adiabatic expansion models are 
very minor, as long as the same asymptotes are reproduced 
( Piranl[l9^ leads to a non-relativistic asymptote different from 
the self-similar asymptote). A more pronounced difference be¬ 
tween possible shell evolutions is that between radiative and 


adiabatic expansion, and between large and small initial mass 
content (the latter discussed separately in section l231 l. In the 
radiative case, blast wave Lorentz factor and radius are still dic¬ 
tated by the energy within the blast wave, only now this en¬ 
ergy is diminishing noticeably due to radiative losses, leading 
to faster deceleration. Calculating the cumulative energy loss 
from standard synchrotron afterglow emission (as covered be¬ 
low in section |3]l for a shell model will typically only add up 
to only a few percent well into the non-relativistic stage, justi¬ 
fying the adiabatic assumption from the beginning of this sec¬ 
tion. However, under certain circumstances, for example when 
the unshocked medium contains many electron-positron pairs 
triggered by promp t emiss ion photons ((Thompson and Madau , 


2 OOOI: Beloborodovl 2002 ). conditions leading to radiative blast 


waves may arise at early times. Observational support for ini¬ 
tially radiative blast waves is offered by Fermi EAT gamma-ray 
detections, w hen interpreted in the afterglow blast wa ve frame¬ 
work (see e.g. Ghisellini et al.l20T^ Nava et al.l2013h . Most of 
the cited shell models incorporate the possibility of a significant 
radiative energy loss term. 


2.2. The non-relativistic Sedov-Taylor-von Neumann self¬ 
similar solution 

A full solution for the blast wave in the radial flow case re¬ 
quires equations for the fluid profile everywhere in the flow, not 
just behind the shock front. For the limiting cases, when T 1 
or yS X 0 , this solution is provided by the aforementioned self¬ 
similar solutions, since dimensional analysis indicates only one 
possible dimensionless combination for the remaining variables 
(i.e. E, Pext, radius r, lab frame time t). These analytical solu¬ 
tions tend to be unwieldy, but sometimes reduce to very simple 
form. An example of the latter is the non-relativistic stellar 
wind (pext r case, where 

P = ^i.rlR)pext, V = r/(2f), p = Pextr^lORh- (4) 


Here r/R plays the role of the self-similarity variable, and we 
have used an ideal gas with poly tropic index 7 = 5/3. Clearly, 
these reduce to the non-relativistic shock-jump conditions when 
r —» and 72 X 1 + as can be seen from eq. |2] The radius R 
is in this case given by = {\27t l[prefR'lef]y^^■ The 


mass dM within a shell at radius r is 167 r(r 2 //?)p„, dr, con¬ 
firming once more that the swept-up mass in the blast wave is 
concentrated at the front. 

The full solution for arbitrary k can be obtained e.g. by gen¬ 
eralizi ng the k - 0 case as described in iLandau and Lifshitz 

(Il959h : 
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Here Cs is the speed of sound. In our case, p - 3pc^l5. The 
self-similarity yariable ^ s r/R is equal to 1 at the shock front, 
and the shock-jump conditions therefore yield y(l) = 3/4, 
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G(l) = 4, Z(l) = 5/16 as boundary conditions for the self- 
similar functions. Plugging the self-similar Ansatz provided 
by eq. |5] into the fluid hydrodynamical equations and solving 
the resulting differential equations, eventually yields (exclud¬ 
ing k - 2, which is a special case where terms drop out of the 
equations early on): 

.5-^ . / V -5 + k + AV y/ 5V-3 Y 

^ iy(i)j i-5 + ^ + 4y(i)j i5y(i)-3j ’ 

5(1 - y)y" 

3(5y-3) ’ 


2.3. The ultra-relativistic Blandford-McKee self-similar solu¬ 
tion 


On the other end of the velocity spectrum sits the 
Blan dford-McKee self-simi l ar so lution for ultra-relativistic 
flow (IBlandford and McKeel Il976h . As already suggested by 
the facts that in this stage f - 1 - 1 / (Ty^) -H 6>(y“"^), and the 
width of the shell AR oc Rfy^, a Taylor-series expansion around 
7 “' i 0 for the fluid profile will typically have y^^, rather than 
y~^ as its first non-constant contributing order. From energy 
conservation within the blast wave that moves at nearly the 
speed of light, one obtains: 
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2(9 - Ak) 


2k 


V6 = -i 


3{k-l)’ k-5 

From dimensional analysis, we obtain for the radius 

2 41 /( 5 -*) 
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where f can be found using energy conservation, leading to 
167r r‘ ,1 9 




.k-5 


(5 - ky 




y 2 + 


10 


Z]^df. 


(8) 


In the ISM case, >6 =» 1.15, in the stellar-wind case, jS « 0.92. 

The late-time non-r elativistic solution pri mary applies to ra¬ 
dio observations (e.g. Waxman et al.l 1998h . since emission in 
higher frequency bands such as optical and X-rays will have 
already dropped below the detection threshold of most instru¬ 
ments. The distance scale at which a typical blast wave be¬ 
comes non-relativistic, is vast (as can be confirmed by plugging 
Py = 1 into equation[3]l. The Sedov-Taylor solution for nonzero 
k provided above, is therefore not likely to occur in nature in 
such a clean fashion. The blast wave has expanded well be¬ 
yond the sphere of influence of its progenitor, making it more 
likely that its current environment is approximately ISM-like, 
or shaped by some com plex interaction between wind bubbles 
from surrounding stars ( Mimica and Gianniosl 2011). Never¬ 
theless, the non-zero k case is relevant as an asymptotic solution 
to long-term evolution of non-I SM hydrodynamical b last wave 
simulations (e.g. those done by De Colle et al.ll20125) . 


17-4k 


(9) 


where F the shock Lorentz factor. The numerical factor in this 
equation follows again from integrating the (rest frame) energy 
density over the self-similar fluid profile (provided below), as in 
eq. El Note that the shocked fluid is relativistic ally hot (as can 
be seen from the shock-jump conditions in eq. |2]), and y - 4/3. 
Since according to eq. |9] F^ oc it also follows that 


R — ct 



2(4 - k)F2 


( 10 ) 


This equation explains the extreme variability of GRBs and 
early afterglow. A photon emitted from that part of the shock 
front that is moving directly towards us, is observed at 


(1 -Hz) “ 2(4 - k)F2 ’ 


( 11 ) 


where z gravitational redshift, and if the observer time is set to 
zero at the point when the GRB is first observed. The Lorentz 
factors c an reach incredibly high values, 100-1000 and beyond 
(see e.g. Racusin et Ii1l2^ . 

Continuing the Blandford-McKee solution, we again com¬ 
bine the shock-jump conditions and hydrodynamical equations 
with the self-similarity Ansatz, this time for self-similarity vari¬ 
able obtain 


= [l+2(4-k)F2](l-r/[cf]), (12) 


where p' expressed in the lab frame. The ultra-relativistic limits 
of eq. |2]are reproduced taking X 1, corresponding to the 
position of the shock front. In contrast to the non-relativistic 
self-similar solution, the relativistic version does not work all 
the way to the origin, but applies only to the relativistic part of 
the outflow, which is where almost all matter and energy reside 
anyway. In numerical simulations, non-physical values can be 
avoided by simply adding 1 to the profile for y^, when setting 
up the Blandford-McKee solution as initial conditions. 


3 







































2.4. Jetted outflow 

A major complication to the simplified picture sketched 
above, is that GRB blast waves are extremely likely to be col¬ 
limated into t wo diametrically oppo site jets with narrow open¬ 
ing angles 6q ( Rhoadsl 1997 . 1999ll . Due to strong relativistic 
beaming, this is not immediately apparent in the light curve, 
as only a small patch dd ~ l/y is visible at first. Unless the 
geometry of the out flow is radically dif ferent from radial flow 
(e.g. cylindrical, see Cheng et aP 2001 ) or the fluid properties 
are strongly dependen t on angle (e.g. as in ‘structured jet’ mod¬ 
els, iRossi et al.ll2002h . it is not possible to tell apart collimated 
and spherical flow at this point. Additionally, if the outflow 
was launched radially, it will take time for causal contact to be 
established across angles and for sideways motion to become 
apparent in the observer frame, postponing a deviation from 
radial flow even for a conic wedge of limited opening angle. 
Sky images are of no help either, since GRBs are typically too 
distant for sp atially resolved observations (the exc eption being 
GRB 030329. lTavlor et al.ll2004 lOren et alj|2004 . 

Jet collimation therefore has to be inferred indirectly, and 
this can be done in various (model-dependent) ways. First, 
one can compare the number of detected bursts to predicted 
rate of occurrence from a given model or to actual detec¬ 
tions of expected count erparts (i.e. supernov ae, in the case 


of long bursts, see e.g. ISoderberg et al.ll2006l) . Second, one 


can compare between early and late inferences of the energet¬ 
ics of the blast waves. At early times, assuming radial flow, 
the relevant energy is the energy per solid angle, or EjsoK^n), 
where iso stands for ‘isotropic equivalent’. If the jet subse¬ 
quently spreads out sideways and becomes spherical (allow¬ 
ing eventually for the non-relativistic radial-flow self-similar 
limit), the energy per solid angle becomes Ejeil{4-n), where 
Ejet ~ EisoO^I2, for small opening angles. Comparing between 
early and late ti me calorimetry should therefore yield openin g 


early ana late ti me calorimetry snouia meretore yieia openin g 
angles (see e.g. Berger et al. 2004 : Shivvers and Berge^boi ih . 


Even if no late time calorimetry is possible, the often extremely 
high values for Ej sn resulting from early-time calorimetry (e.g. 


Cenko et al. 2010ll . already hint that the actual energies in the 


jets are likely lower. When collimation is accounted for (and 
jets with typical jet half opening angles of 6'’ are inferred, see 
e.g. Racusin et al, 20091: Rvan et al. 201^ the re sults tend to 


cluster around Ej^, ~ 10^' erg ( Bloom et al. . 2003h . 


The third way of inferring collimation, is by looking for sig¬ 
natures directly in the light curve. Two effects will lead to a 
steepening of the temporal evolution. On the one hand, the vis¬ 
ible patch will grow as relativistic beaming weakens and re¬ 
veal a lack of emission from beyond the edges once y l/0o- 
This effect strongly depends on viewing angle, which puts the 
visible patch init i ally closer to the jet edge if moved off-axi s 


( van Eerten et al.[ 2Q10a : van Eerten and MacEadvent 2012b ). 


On the other hand, the jet will start to spread out sideways. The 
over-pressured edges (relative to the circumburst environment) 
will do so immediately, and as the fact of the emergence of 
the jet into its environment is communicated towards the jet 
axis, more of the jet will follow. The broader jet will sweep 
up more material per unit time, leading to a faster decelera¬ 
tion. Because the spreading velocity of the jet is suppressed by 


a factor 1/y in the observer frame, this sets off a feedback loop 
where a slower jet is seen to sweep up even more material. Due 
to this 1/y suppression factor from the Lorentz transform be¬ 
tween the frame comoving with the blast wave and our frame, 
the sideways spreading of the jet becomes noticeable again once 
y ~ 1 / 00 - If the jet were in fu ll causal conta ct, the spreading 
behavior would be exponential (lRhoadsLll999l) . 

No exact analytical solutions exist for spreading jets, even 
when starting from a ‘top-hat’ conic wedge out of the spher¬ 
ical self-similar solution (once sphericity is dropped, ini¬ 
tially structured and multi-component jets can be ass umed 


too, see e.g. Berger et al. 2003 : Rossi et ^ 2002), al 


Rhoadsl 19991 Sarietal. 1999: Kumar and Panaitescu 

2000 

Huang et al.l2000' van Eerten et al.l2010a: 

Wvgoda et al. 

2011 

Granot and PiranI 2012h. Unfortunatelv, 

approximate mod- 


els tend to be notoriously sensitive to the choice s made for 
the underlying simplifications (see iGranoti l2007l for a dis¬ 
cussion). The picture that has emerged from relativistic hy 


drodynamics (RHP) sim u lations ( ^ang and MacEadvenl2009 


van Eer 

ten et al. 201 Oat 

Wvgoda e 

ral.1 2011: 

De Colle et al. 

201 2 hl 

van Eerten and MacEadvenI 

20131). is one where, for 


realistic opening angles (where 6q 0.05 rad does not 
apply), jet spreading does not achieve the runaway behav¬ 
ior and exponential increase in opening angle expected in 
the ultra-relativistic limit, but sta ys closer to logarithmic 
dvan Eerten and MacPadveiil 2012bll . A key reason for this is 
that full causal contact along all angles of the fluid takes time to 
establish, leaving little time in practice for an ultra-relativistic 
spreading regime because the blast wave quick ly becomes 
trans- relativistic following the onset of spreading (Ivan Eerten . 
20131) . As a result, the effect on the light curve following the 
‘jet break’, is due to the joint impact of both the edges be¬ 
coming visible and the onset of spreading, with neither over¬ 
whelming the other. Post-break slopes remain dependent on 
observ er angle and can b e used as means to constrain jet orien¬ 
tation ( Rvan et al. . 2015 ). 


2.5. Reverse shocks and injection of energy into the flow 
The preceding sections describe the subsequent evolution of 
an instantaneous explosion without initial mass content, i.e. a 
single forward shock moving into the circumburst medium. But 
unl ess the jet is dri ven nearly completely by pointing flux (see 
e.g. Lvutikovl2006l) . a certain amount of initial mass is expelled 
with the explosion. A simple way of incorporating this would 
be to add a mass Mq to the shell model (but prior to deriving eq. 
[2 because the ejected mass is presumed to reside in a cold shell 
where all energy is converted into kinetic energy, while eq. |2 
only connects mass to shock-jump conditions). The presence of 
initial ejecta mass will postpone deceleration of the blast wave, 
which will initially coast along at fixed velocity in ballistic mo¬ 
tion. A forward-reverse shock system is established, with the 
forward shock moving into the circumburst medium and the re¬ 
verse shock heating up the ejecta. The width of the ejected shell 
influences the dynamics as well. Even an initially infinitesi¬ 
mally thin shell will stretch out to width AR ~ Rjf^, similar to 
a decelerating blast wave. Only for shells wider than this, the 
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initial shell width has to be taken into account explicitly when 
computing the deceleration radius, marking the turning point 
between coasting and decelerating of the blast wave. The two 
types of shell, wider or narrower than the intrinsic blast wave 
width, have bee n labeled ‘thick’ and ‘thin’ shells respectively 
in the literature ( Sari and Piranl 1995h . For thick shells, the re¬ 
verse shock achieves relativistic velocity (in the ejecta frame) 
before crossing the ejecta. The pre-deceleration stage and re¬ 
verse shock crossing were expected to be visible briefly (mainly 
in the optical), during the early evolution of the afterglow, for 
about 10^ s. 


However, the Swift satellite (iGehrels et al.L 120041) has re¬ 
vealed instead the existence of an extended (10^^ s) early, flat 
phase of X-ray d ecay to be a common occurrence in a fter- 
glow light curves ( Nousek et al. . 2006t Zhang et all 2006 ). A 
similar sta ge has been found to exist in the optical as well 
(see e.g. Panaitescu and Vestrand 2011 : Li et al. 20121) . al¬ 
though a joint explanation for both X-rays and optical is com¬ 
plicated by the existence of very early steep decay in the X- 
rays (likely connected to the end of the prompt emission, rather 
than the afterglow), complex temporal optical emission profiles 
(maybe multi-component emission, from both reverse and for¬ 
ward shock) and the close proximity of the end of the shal¬ 
low de cay phase and of what is typically interpreted as a jet 


break (iLi et al 


gle effects (e.g^ 


2015h. If not explained from viewing an- 


Eichler and Granot 20061) or evolvin g micro¬ 


physics (e.g. iGranot et al.l l2006t iFilgas et al.l 1201 Ih . plateau 


stages point to some form of prolonged energy injection into 
the ejecta, where a continuous ‘push’ from the back delays de¬ 
celeration. This can take various forms, such as late catching 


system (e.g. 

Panaitescu et al.l 1998: Rees and MeszarosI 1998: 

Sari and MeszarosI 

I2OOOI). 

long-term source luminosity (e.g. 
or conyersion of Poynting flux from 

Zhans and Meszarosl2001 

the eiecta fe.g. Usoyll992l 

Thomosonll994l). In afterglow anal- 


ysis, this injection can be mo deled in the form of a power law 
increa se in ejecta energy (e.g. Zhang et^l2006t Racusin et al 


2009h . In the case of a relativistic reverse shock (in the frame of 
the inflowing material) and gradual, po wer-law type injection , 
one can even ma i ntain se lf-similarity (IBlandford and McKee . 
19761: van Eerten , 2014ah . Another promising modeling ap¬ 
proach which has been applied to data directly, is dropping self- 
simila rity, but ma i ntaining a simplified descriptio n for the late 
shells ( Uhm et al. . 2012t De Pasauale et al.l 2015). It should be 
emphasized that, although long-term engine activity is certainly 
a possible explanation for these early st ages (potent i ally requir¬ 
ing a magnetar engine rnodel, see e.g. Usov|ll992 ; I Thompson 
1994t Dai and Lu 1998 ; Zhang and Meszar^ 200 ih. iet break¬ 


out is a messy even t (see e.g. IWaxman and MeszarosI 12003 


Morsonv et al. 2007h which might well naturally introduce an 
extended observable early stage even for briefly active engines, 
before moving toward s the asymptotic regime of a de celerating 
relativistic blast wave ( Duffell and MacEadvenl 2014). 

A potential means of distinguishing between engine models 
are the correlations fou nd between plateau e nd times and lu- 


min osity in the X-ravs (IDainotti et a ll 12008) , and in the opti¬ 


cal dPanaitescu and Vestrandl 1201 It iLi et al 


20121) . The two 



Figure 1: Typical synchrotron spectrum in the slow cooling case. 


correlations have different slope, which is consistent with the 
optical and X-ray emissio n typically be i ng ob served to be in 
different spectral regimes (IGreiner et all 1201 ll) . They emerge 
naturally from a synchrotron forward-r everse shock system 


( Leventis et al. . 2014 : van Eerten . 2014ah. but require the pres- 


ence of a relativistic reverse shock dvan Eerten! 2014bl) . or, 
‘thick’ shells rather than ‘thin’ shells. It is not clear how strong 
the emission from the reverse shock will be in reality, since the 
relative strength of reverse shock emission is sensitive to model 
assumpti ons, such as the degr ee of magnetization of the ejecta 
(see e.g. Mimica et al. 2009bl) . a nd can vary wildly e ven for a 
standard synchrotron model (e.g. Leventis et al.ll2014 ). 


2.6. Further complications 

As already alluded to above, the initial geometry of the ejecta 
and the environment of the burster provide two obvious compli¬ 
cations to the standard picture. Even the direct environment of 
the burster can reasonably be expected to have a complex shape. 
The stellar wind profile will only extend to a finite range and be 
influenc ed by photo-ionization, stellar rotation and fluid insta 


bilitie s (|Eldridgg-gLal 


2007: yan Marie et al 


2006; van^arlee^al^ 2006; Eldridge . 


2008; yan Marie and Keppensl 20121) . 


Late time mass loss of the progenitor system is likely erratic 
(IMesler et al.L l2012h . Although circumburst mass transitions 
are not expected to introduce sudden changes in the obseryed 
light curyes from the forward shock (Nakar and Granot . 20071 : 
yan Eerten et all 20091: Gat et all 2013 ; Geng et alT 2014 ), an 


oyerall slope transition can reasonably be expected, which 
could explain k me asurements other than k - 0 or k - 2 
(ICurran et al.L l2009h . Additional emission might be generated 
by a complex shock structure following m ultiple interactions 
( Uhm and Zhang , 2014 : Mesler et al.l 2014 ). 


3. Emission 

The synchrotron spectrum consists of a series of connected 
power laws, separated by break frequencies, and eyolyes 
in a characteristic manner during the lifetime of the blast 
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wave ( Meszaros and Re^ 1997 : Wiiers et al. . 1997 : Sari et al 

1998h . An example is provided by Fig. [T] We will discuss this 
slow cooling case first. A shock-accelerated electron popula¬ 
tion neije) with power-law p ~ 2.3 is typically assumed, with 
neije) and je expressed in the comoving fluid frame. 

The lower cut-off value for this distribution is If we further 
parametrize the shock-accelerated electron spectrum using 
(typically taken ~ 1), the fraction of electrons that get shock- 
accelerated, and Ce (typically ~ 0.1), the fraction of available in¬ 
ternal energy in the fluid that goes into the non-thermal electron 
population, one can derive (by equating integrals over electron 
number density distribution and electron energy density distri¬ 
bution ye«c(Tc)'«cC^ to available total number density and en¬ 
ergy density respectively): 


Tm 


2 — p I eg e tnp 


(13) 


1 - p \^N pnieC^ 

where lUp proton mass and electron mass. If the power law 
distribution slope p is too shallow {p < 2), one can choose to 
either maintain a physically plausible proportionality between 
jm and ejn, or have the upper cut-off to the particle popula¬ 
tion (which can be ignored for p > 2) d ictate instead, in 


order to maintain the interpretation of (iBhattacharvai 12001 
Dai and Chengl 200 ih . 


According to synchrotron theory, the local spectrum from an 
electron at energy je peaks around 


3 

T^e -- 

47t nigC 


(14) 


in the frame comoving with the fluicd and where qg electron 
charge. Magnetic held B is typically parametrized again via a 
fraction of available energy, according to l{^n) = egc (and 
with cb typically ~ 0.01). The critical frequency v,„ shown in 
Fig. m represents the average critical frequency for the com¬ 
bined emission of all local synchrotron spectra and their local 
values, and expressed in the observer frame. The spectral 
slopes of 1 /3 and (1 - p)l2 at both sides of also follow from 
stand ard synchrotron theory (see e.g. 

1979h . 


Rvbicki and Lightman 


The dependency of the Hux on the model parameters Eiso, 
Pext, SB, Se, z, di (luminosity distance) in a given spectral 
regime can now be determined as follows. The emission coefh- 
cient peaks according to synchrotron theory at 


1 ^ql 


^base 


nieC^ 


^NnB, 


(15) 


in the frame comoving with the huid and where n the local co¬ 
moving huid number density (such that ^ffn the number density 
of non-thermal electrons). On both sides of v'„, we have 




fv = 


^base 


^base 


1/3 


— , V < V™ < V. 


(i-rt/2 


V„ < V <Vg 


(16) 


^This choice of notation was made for consistency with the literature. Note 
that we now have v' and p in the frame comoving with the fluid, while v and p' 
are expressed in the lab frame. 


where V the observer frequency and v' the cooling break in the 
frame comoving with the fluid, which we will discuss below. In 
our frame y^e'/, as the dependency evolution of flux on 
model parameters will be dictated by emission from material 
moving (nearly) straight towards the observer. The observed 
flux is then 

Fv',@ oc (17) 

“l 


Here V is the emitting volume, the product of area (R/y)^ 
and depth AR oc Rjj^. The Lorentz factor in the area re¬ 
flects the size of the visible patch due to beaming. Without 
sideways spreading, it would be sufflcient to omit this fac¬ 
tor in order to obtain post jet-break flux values. Spreading 
models quickly become more complicated (simulation-based 
po st-break light curve slopes for th e ISM case are provided 


by Ivan Eerten and MacFadvenll2013L These simulations reveal 


steep temporal slopes of ~ -2.7 above v^, once the transi¬ 
tion to the post-break regime has completed, which is strik¬ 
ingly steeper than indicated by the Swift sample). Observed 
frequencies are related to comoving frequencies via the usual 
Vg » yv'/(l + z), and observed time and emission time via 
eq. [TT] Using the jump-condition values from eq. |2] and the 
dynamics from eq. |3l is then sufflcient to determine tfle ex¬ 
act dependence of flux on the model parameters. Synchrotron 
flux equations for afterglow blast waves can be found at vari¬ 
ous places in the literature, including extensions such as trans- 


—-,- 

e.s. iMeszaros and Rees 

19971 

Sari et al. 

1998; 

Waxman et al. 

19981: Gruzinov and Waxman 

Il999l Granot and Saril 12002; 

van Eerten and Wiiers 

20091 

Leventis et al. 2012: Yi et al. 

20131: Gao et al.ll2013: van Eertenll2014ah. 


If not all available electrons are shock-accelerated into a 
non-thermal population, i.e. < 1, some will remain in a 
Maxwellian distribution. This can potentially impact the syn¬ 
chrotron light curve. Unfortunately, can not be derived from 
observations of the power-law electrons directly, even if all 
spectral regimes could be observed. A full degeneracy between 
and other model parameters exist, where a decrease in 
can be compensated for by a simultaneous linearly proportional 
decr ease in e,,, cb and linearly pr oportional increase in Eiso and 


Pext 


( Eichler and Waxmanl 20051) . 


3.1. Electron Cooling 

Another characteristic of the synchrotron spectrum is pro¬ 
vided by electron cooling, since energetic particles use their 
energy very quickly through synchrotron emission: 


dye _ 4o-7-y^ jg An 

df' hnigC 3n dt' ’ 


(18) 


where ctj- the Thomson cross-section and f in the comoving 
fluid frame. Beyond the cooling break Vg, the effect of cooling 
becomes important and the first term on the right above (the 
synchrotron loss term) dominates. Below the cooling break, 
the other term (cooling through adiabatic expansion, note the 
‘stretching’ h/n) dominates. The cooling break Vg can lie either 
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Figure 2: A fast cooling synchrotron spectmm. 


above or below v^, yielding spectra designated as slow- and 
fast-cooling respectively. 

A simple approaching to modeling the behavior of the cool¬ 
ing break, is to ignore the spatial structure of the blast wave 
for the purpose of determining the cooling time, using a steady- 
state approximation to eq. [T8]and equating overall cooling time 
to burst duration, resulting in 


6nmeCy 


(19) 


connected to a characteris t ic freq uency in the usual manner 
via Vc oc dSari et al.L 1199811. This approach has been 
used in simula t ions as well (e.g. Izhang and MacFadvenlE009 


De Colle et al. 2012bt van Eerten et al. 2012h . even though 
other quantities (vm, peak flux) are calculated completely lo¬ 
cally. The correct scalings and temporal evolution are repro¬ 
duced in this manner, but this hybrid approach (in the simu¬ 
lation case), implies an offset of the cooling break relative to 
a fully local approach to cooling. The reason for this is that 
the full fluid profile provides a dimensionless constant of in¬ 
tegration when computing the cooling break from local spec¬ 
tra (which exhibit an exponential cut-off locally, rather than a 
power-law transition, only globally adding up to such a tran¬ 
sition) that is different from the one provided by an effec¬ 
tively flat fluid profile. The good news is that this off-set 
remains ess entially constant throug hout the evolution of the 
blast wave ( van Eerten et al.l 2010ali . but the effect on multi 


band analysis of afterglows can be substantial (iGuidorzi et al 


2OI4I) . A local approach to electron cooling in simulations 
(e.g. by rewriting eq. [18] in the form of an advection equa¬ 
tion for Ym), requires extre me resolutions, which is c hallenging 
already in one dimension (Ivan Eerten et al.L 1201 Obh . but only 
achievable in multi dimensions by specialized methods (e.g. 
van Eerten and MacFadvenll2013 ). 

In the global approach, the emission coefficient equations [T6l 
can be extended to include the effect of electron cooling as in- 




mpeak 



Figure 3: self-absorbed spectrum beyond v,„, in terms of regular ordering for 
Va < Vm- 


dicated in Figs. [T]and|2| leading to 


^base 


^base 


1/2 


-p/2 


<v'< V, 


V, < y < v„ 


( 20 ) 


Here tracks the peak of the spectrum, which not neces¬ 
sarily coincides with v,„. In the local approach to cooling, this 
only describes the shape of the globally emergent spectrum (i.e. 
flux, rather than emission coefficients). The -1/2 power can be 
understood as follows. The frequencies in this regime probe 
electron Lorentz factors below the injection value of Ym, while 
the cooling timescales are extremely short since we are above 
Vc (meaning that the shape of the injected profile is not relevant, 
so no p-dependency). All available energy is quickly radiated 
away. Per frequency, this yields Fy oc Yei'fieC^/[ yIB] oc oc 
v- 1/2 


3.2. Synchrotron Self-absorption 

At low frequencies (typically in the radio regime), syn¬ 
chrotron self-absorption (ssa) becomes important and a charac¬ 
teristic break Va appears in the spectrum, below which the flux 
drops off more steeply. Since ssa is not a scattering process, it is 
relatively s traightforward to model in a linear radiative transfer 
approach, ( Granot et al. . 1999[ Granot and Saril 2002 ) that can 
also be applied to adiabatic relativistic blast w ave simulations 
( Mimica et al.l 2009a; van Eerten et al.L 2010b ). In the absence 
of electron cooling, the absorption coefficient due to ssa is given 
by 


a'y’ ~{p-l)(p + 2)fivn 


(vr^i^ 

y„,167rm2c2 \v'^ 


( 21 ) 


where k - 1/3 if v' < and k = -p/2 otherwise. In a simpli¬ 
fied computation in order to obtain the model parameter depen¬ 
dencies, we need to consider an emitting surface A cx (RlyY 
and a source function below Vg, rather than the emitting vol¬ 
ume f or the optically thin case ( Sari et al. , 1998 ; Waxman et all 
19981) . i.e. 

Fy,,^oc^A — . 

Cl I Cly' 


( 22 ) 
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In the lab frame ay> =» o'ylj- In the slow cooling case and with 
Vm > Va, the spectral slope 2 below Va (see Fig. [T]l follows 
from a comparison between eqs. |2T] and [T^ In case and 
Va flip, a new slope of 5/2 is introduced, as can be seen from 
the same equations. Once Va > v,„, the temporal evolution of Va 
will change, as will the peak of the spectrum, which now occurs 
at Va, rather than v^. Nevertheless, these differing characteris¬ 
tics follow from the regular self-absorption and evolution, 
as can be seen geometrically in Fig. |3] where Va,i the position 
when extrapolating the Vm > Va case and Vg s the actual posi- 
tion of the self-absorption break (see also iLeventis et al.ll2012 
for a discussion. The numbers 1 and 5 were chosen to match 
the notation from Granot and Sarill20()2 ). The model parameter 
dependency and evolution of Va (i.e. Vaj) can be determined by 
looking for the meeting point between the flux according to eq. 
[22] and eq. [17] using v < Vm for both. 

As said, the case Va,$ > Vm does not introduce anything new, 
although the flux equations look different. This minor observa¬ 
tion has the practical implication that, when constructing scale- 
invariant spectral templates from high-resolution simulations, 
all that is needed are the temporal evolution curves for v,„, Vc, 
Vflj and Fpeak, with the understanding that the latter does not 
coincide with the actual spectral peak once s > v,„. 

The fact that self-absorption renders only the front of the 
blast wave visible to the observer should serve as caution when 
attempting to seek out early time emission from a reverse shock 
in the radio domain: depending on density and profile of the en¬ 
vironment, this component might well be hidden from view by 
ssa. A precise analysis of the early stage radio emission is fur¬ 
ther complicated by the role of electron cooling. At early time, 
we might also be observing the fast cooling case, rather than the 
slow cooling case, and the cooling break will help shape the ab¬ 
sorption coefficient. At this point the difference between local 
and global cooling emerges again as well. An exact treatment of 
local cooling will actually introduce additional spectral regimes 
that are not apparent in a global appr o ximat i on. These topics 
are dis cussed further by Granot et al. ( 1999h : Granot and Sari 

(I 2 nn 2 h . 


3.3. The sharpness of spectral breaks 

The connected power law description of the synchrotron 
spectrum that we applied above, is of course an approxima¬ 
tion. In reality, the asymptotic regimes approach one another 
smoothly, with the underlying shape for a single electron spec¬ 
trum being an integrated modified Bessel function of fractional 
order. Since the full expression for a spectral transition is cum¬ 
bersome and, more importantly, since the spectral sharpness 
also depends on the fluid configuration that shapes all the local 
contributions to the emergent spectra, it is more convenient in 
practice to use approximate formulae, typically smooth power 
laws of the type 


Y(x) = To 




(23) 


and varieties thereof, with the most useful shape depend¬ 
ing on e.g. the sign of the transition 02 - ai. The larger 


s, the sharper the transition. Such approximations can be 
applied to temporal transitions (e.g. jet breaks) as well and 


Beuermann et al. 

1999; 

Harrison et al. 

19991: 

Granot and Sari 

2 OO 2 I: van Eerten and Wiiers 2009 

Leventis et al.l 2012: 


van J 


Eerten and MacFadven 2013li . Smooth power laws can 


lead to significantly better fits to data, and the flux exactly at 
a spectral transition can differ up to an order of magnitude 
from sharply connected power laws. The exact value of s is 
hard to determine in practice from the data, and good fits can 
typically be obtained for a range of values. The theoretical 
values of s are influenced by many things. In addition to 
the fluid profile, the closeness of other spectral transitions 
also plays a role even if their presence is not immediately 
apparent from the data. Finally, the sharpness of the spectrum 
is sensitive to the orientation of the jet, since a given observer 
time c orresponds to a different set of em ission times for each 
angle (Ivan Eerten and MacFadverul2012ah . 


3.4. Further complications 

Aside from the complexities introduced by more realistic 
particle spectra and by multiple emission sites for synchrotron 
emission (e.g. a forward shock and a reverse shock), the 
most obvious further complication is the addition of other ra¬ 
diative processes. Of these, inverse Compton scattering and 
synchrotron self-Compton scattering are the most likely can¬ 
didates. Inverse Compton scattering has a noticeable i mpact 
on the cooling of elec t rons when cr Ce ( Sari et al.l 1 19961 : 
Panaitescu and Kumail 200(i Sari and Esinl 12001 ). This ef¬ 
fect can be included in prescriptions for synchrotron spectra. 
Even when not observed directly, Inverse Compton scattering 
will shift the cooling break downwards, a nd, in the fast cooling 


stage, the self-absorpti on break upwards. (ISari and Esinl 12001 
Granot and SariL 2002). 


Other factors that complicate interpreting the observed emis¬ 
sion in terms of a synchrotron blast wave origin, is the con¬ 
tribution from completely separate components, such as host 
galaxies and supern ovae, or dust echoes of the prompt emission 
( Evans et al. . 20141) . 


4. Model-based data analysis 

When analyzing the data, a number of approaches can be 
taken. The conventional approach has been to start from anal¬ 
ysis of light curves and spectra (when available) in terms of 
simplified heuristic fit functions, typically power laws, aug¬ 
mented where appropriate by descriptions of extinction and ab¬ 
sorption due to intervening material and host galaxy and su¬ 
pernova flux contributions. This approach results in a con¬ 
cise description of the data (in itself already useful) that can 
subsequently be interpreted in terms of physical models, un¬ 
der the assumption that the essence of these models can be 
captured sufficiently in the form of power laws. On the other 
hand, once can test physical models against the data directly 
in a manner tha t does not require synchronous mul t i-band ob- 
servations (e.g. Panaitescu and Kumar 2001 , 2002; Yost et all 
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Figure 4: Evolution of the synchrotron spectral peak, computed from lah- 
frame (blue) and boosted-frame (red) simulations. Dashed grey lines denote 
the asymptotic self-similar limits. The drop between 10® - 10* days describes 
the effect of a jet-break. The additional bump at late time the rise of the coun¬ 
terjet. 


299?) 

which has become increasingly popular in recent years 

(e.g. 

van Eerten et al. 2012: Leventis et al. 20131: Laskar et al. 

20131 

Perlev et al.l2014 ; Urata et al. 

20 i4 

Guidorzi et al.l2014: 

Rvan et al. 20151: IZhane et al. 2014 

). Rather than yielding the 


‘best’ short-hand description of the data set, this immediately 
addresses the question of whether a given physical model can 
explain it, and to the extent that it can, provides estimates for the 
model parameters. An additional advantage of this approach is 
that arbitrarily complex light curve shapes (as provided by the 
underlying model) can be accounted for. 

These model fitting approaches also naturally connect to 
probabilistic data analysis methods, including Bayesian in¬ 
ference, which are having a transformative effect on the 
field. Software packages that implement numerical meth¬ 
ods such as affine invariant Mark o v-Chain Monte Carlo 
(MCM C, Goodman and We are 201(i ForemamMactevet^ 


2013h or mu l ti-mod al nested sampling Jperoz et all 


2009; 


Buchner et al.L 1201 4) are becoming readily available and are 


ideally suited to GRB afterglow analysis, due to their capabil¬ 
ity of dealing efficiently with expensive fitting functions, large 
numbers of free parameters and bimodal posterior distributions. 
The possibility in Bayesian data analysis of marginalizing over 
nuisance parameters is extremely valuable for assessing the per¬ 
formance of a model in case not all parameters can be con¬ 
strained fully (as is often the case for afterglows, since this 
would require broadband data covering all spectral regimes of 
the synchrotron spectrum; such data sets exist, but are still rare, 
although the sample is growing). 


4.1. Simulations 

Because the deceleration to trans-relativistic veloci¬ 
ties and the spreading behavior of relativistic blast waves 
are difficult to capture analytically, many groups have 
used RHD simulations in one and more dimensions to 


study this crucial stage of aftergl ow blast wave evolu¬ 


tion (including, but not limited to Kobavashi et al.l 119991: 


2007 


Zhane and MacFadven 

2011 


Wvsodaetal. 2011 


200 ^ 


van Eerten et al 


De Colle et al 


2012 a: 


van Eerten and MacFadven 2013: Duffell and MacFadvenI 


2013h . Most multi-dimensional simulations employ special 


strategies to deal with the stringent resolution requirements (i.e. 
AR ~ R/F^ with F » 1), specifically adaptive-mesh refine¬ 
ment, where the grid resolution is locally dynamically doubled 
(or halved) along each dimension, depending on the variability 
of the flow. With this approach, the six orders of magnitude be¬ 
tween initial shell width and final shell radius could be more or 
less covered numerically. However, it is important to note that 
all simulations prior to 2013 failed to fully resolve numerically 
the blast wave at the key dynamical stage of jet spreading. That 
qualitatively correct behavior w as nevertheless reproduced, 
could only be confirmed recently (Ivan Eerten and MacFadven , 
2013h by simulations in a specialized Lorentz-boosted frame (a 
natural antidote against 1/F^ related issues). These resolution 
issues also impact light curves computed from simulations. The 
strategy for such computations is not dissimilar from inferring 
flux from simplified models: take the dynamics as a given, 
assuming adiabatic expansion, and employ radiative transfer 
or, above the self-absorption break, perform a straightforward 
summation of locally emitted power. As described previously, 
the time-dependent synchrotron spectrum can be captured 
concisely based on its key characteristics. An example is 
provided by Fig. |4l showing both lab-frame and boosted frame 
evolution curves for the spectral peak. 

Because derived synthetic light curves contain emission from 
many different emission times arriving at each single observer 
time, the resolution issue is in practice not problematic at the 
light curve level for homogeneous medium simulations: the 
resolution issue is at its most severe early on during the simula¬ 
tion, and the emission from these times is observed jointly with 
emission preceding the simulation starting time (i.e. from the 
analytical self-similar Blandford-McKee solution that provides 
the initial conditions for the simulation). Only when jets be¬ 
come too narrow (0o 0.05 rad), observer frequencies drop too 

far below at early times, or for certain circumburst medium 
profiles (including, unfortunately, stellar wind), this becomes a 
potential issue. 


4.2. Scale-invariance and model templates 

The complete evolution of a blast wave, from early time 
conic wedge out of the relativistic self-similar solution to the 
spherical non-relativistic self-similar solution, exhibits a use¬ 
ful scale inva riance that can be emp loyed for model compar¬ 
isons to data ( van Eerten et all 20121) . This invariance goes be¬ 
yond self-similarity and follows straightforwardly from dimen¬ 
sional analysis. The fluid equations can be expressed in terms 
of spacetime coordinates = rcjt, f, 6, rather than r, f, 9. In 
the non-relativistic limit i 0, while 9 no longer applies due to 
sphericity, and the self-similar solution of single variable ^ is re¬ 
covered. In the ultra-relativistic limit 1, while 9 again does 
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not apply due to the radial flow assumption, leaving again a self¬ 
similar solution. Note that x - x(^) the Blandford-McKee 
self-similarity variable is more practical in this limit. But even 
for intermediate values of it remains true that rescaling in 
explosion energy or circumburst density can be compensated 
for with a rescaling of the coordinates. Taking Eiso kEiso, 
PrefRlj ^ AiprefRlf), ^ ^ and t ^ 

leads to the invariance ^ ^ ^ and 6 6: a big¬ 

ger explosion explosion (or one in a more dilute medium) goes 
through the exact same stages as a smaller one, albeit at larger 
radii and at later times. In terms of dimensions, we scaled grams 
by a factor of k, and centimeters and seconds both by a fac¬ 
tor {k !The implied scalings for mass densities, en¬ 
ergy densities and pressure are identical: p ^ ^-klO-k)yio-k) 
e —> K^k/{3-k)^3/o-k)^^ p i^-kio-k)^io-k)p_ Lorentz factors 
remain unaffected, 7 ^ 7 . 

The big practical relevance of this comes when building a 
template set out of simulations, for comparison to observational 
data. Two dimensions in parameter space are now accounted 
for (£,50 and p^xt, leaving only Oq). And, although we have no 
exact analytical solution for the spreading stage, the blast wave 
nevertheless segues smoothly from one asymptote towards the 
other. This means that we can use our intuition from the self¬ 
similar asymptotes (with features such as e.g. A£ oc £/r^) to 
guide us towards a suitable compression algorithm for simula¬ 
tion data, rendering the construction of templates feasible even 
if the original simulations are very computationally intensive. 
In this manner, model results along the remai ning Oc dimension 


in par ameter space can be tabulated concisely (Ivan Eerten et al. 


2012h . It should be noted that introducing energy injection does 
not break dynamical scale invariance, even though it introduces 
extra dimensions in para meter space, such a s energy injection 
duration, that scale along ( van Eert^ 2014bh . 

Although dynamical templates are in principle sufficient for 
simulation-based afterglow analysis (with the caveat that radia¬ 
tive transfer based on dynamical templates needs to be cal¬ 
culated on-the-fly), we can do better if the radiative process 
of interest is sufficiently simple: a convenient feature of the 
power law nature of the synchrotron spectrum is that it renders 
it scale invariant ( van Eerten and MacFadvenl 2012a ). Like the 
dynamical scale-invariance, this follows directly from dimen¬ 
sional analysis, and although in synchrotron spectra a number 
of additional constants appear (such as nie), this invariance is 
not compromised within a single power law regime. In the flux 
equation, the role of fractions es, Q and does not change over 
time either. All this implies that evolution curves for and 
other characteristics of the synchrotron spectrum (peak flux, 
cooling break, self-absorption break), even when derived from 
high-resolution multi-dimensional simulations, can be rescaled 
between model parameters, allowing for a synchrotron spectral 
template-based approach to afterglow fitting. 

4.3. BOXFiT and scalefit 

The scale-invariances described above have been used 
to prepare synchrotron templates directly from r el ativis 


tic hydrodynamics simu l ations (^ yan Eerten et al.L 12012 


yan J 


Eerten and MacFadyenl 2012a : [Ryan et ^ , 2015 ). A 


071020 



Figu re 5: A spec t ral tem plate-based fit to the X-ray afterglow for GRB 071020, 
from iRvan et al.l J2015h . Top figure shows the best fit, bottom figure the pos¬ 
terior probability distribution projections for jet opening angle and observer 
angle, marginalized over the other fit parameters (£'/so, pext, ^e)- In a 

single band fit, most parameters remain unconstrained. 
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simulation-based analysis code utilizing a set of 19 homoge¬ 
neous medium simulations an d dynamical scale invariance 
(boxfit, Ivan Eerten et alJ 2012h . is freely available for down- 
loac0. Its follow-up, using scale-invariance at the level of 
spectral templates directly, is currently in preparatio n (scalefit; 


a first application to Swift XRT data can be found in iRvan et al 


2015[ limited to homogeneous medium simulations) and 


will be available in the near future via the same website 
as BOXFIT. The full release of scalefit will include spec¬ 
tral templates based on approximately 70 -f high-resolution 
simula tions including ‘boosted Lorentz f rame’ based simula¬ 
tions ( van Eerten and MacFadvenl 2013 ) and a stellar-wind 
environment. 

When emission is computed from radiative transfer through 
a fluid dynamical solution provided either analytically, from 
simulations directly, or from templates, the resulting spectra 
will be smooth automatically. Direct reconstruction of spec¬ 
tra from scale-invariant templates, on the other hand, requires 
that this smoothness is either ignored or accounted for explic¬ 
itly: spectral sharpness itself is not scale-invariant, because it 
connects two regimes where the flux equations scale differ¬ 
ently for changes in Eiso and pe^t- Additionally, there is the 
practical issue of how to deal with crossings of spectral breaks 
( Leventis et al. . 20121 ) . 

Another complication is that the physical models describe 
spectral flux, while detectors such as Swift/XRT count photons. 
A transition from c ount space to flux space needs to be made 
( Evans et al. . 2009h . while at the same time absorption (in X- 
rays) and extinction (in the optical), need to be accounted for. 
This transition is not completely trivial, and for example re¬ 
quires that some choice be made for the underlying spectral 
shape in the relevant spectral range. Ideally, the spectral shape 
is provided by the physical model as well, and the fit is essen¬ 
tially done in count space. The disadvantage of this approach is 
that this renders the fit software instrument-specific. 


and more complicated initial jet profiles will (at least initially), 
give rise to more complicated light curves. 

The two codes mentioned by name, boxfit and scalefit, of¬ 
fer the means to perform this type of simulation-based model 
fitting, but much work remains to be done. An extensive set 
of templates is being generated including stellar-wind type cir- 
cumburst profile, and the public release for scalefit is in prepa¬ 
ration. The strengths of these model-based codes are that 
they allow for direct model testing and fitting of complicated 
light curve shapes, incorporate advanced statistical methods 
and avoid the need for simultaneous broadband observations. 
The obvious disadvantage of direct model fits, is that features 
not included in the model are not tested for, and might drive 
the fit into a wrong region of parameter space when occurring 
in the light curve and when the model is forced to account for 
them (e.g. interpreting a plateau stage as regular decay with low 
p value). The open source nature of the software should help 
facilitate the addition of additional dynamical models and ra¬ 
diative processes. This can be done both on the level of the the 
dynamical templates (e.g. by including a more complex evo¬ 
lution of the shock-accelerated electrons, Sironi and Giannioi 
2013I) . and on the level o f the spectral templat es (e.g. shifting 


the cooling break directly, Guidorzi et al.ll20l'4 ). 
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5. Conclusions 

Gamma-ray burst afterglows provide a unique opportunity to 
study highly relativistic flows and fundamental plasma kinetic 
theory in a setting that is impossible to repeat on earth. Their 
core ingredients, relativistic blast waves and non-thermal emis¬ 
sion from shock-accelerated particles, can be well understood 
from simplified models, a few of which have been described 
in detail in this paper. Nevertheless, GRB afterglows exhibit a 
rich range of features that follow both from the interplay of the 
standard model components and from the inevitable complica¬ 
tions that can be added to the basic models. Jet and blast wave 
dynamics present an example of this, and it has only recently 
become possible to accurately model the evolution of GRB out¬ 
flows, through high-resolution numerical hydrodynamics simu¬ 
lations, in a manner that can be compared directly to broad¬ 
band observations. Even then, ‘accurate’ has to be understood 
in terms of the simplicity of the initial assumptions for the jets. 


'^HTTP://cOSMO.NYU.EDu/AFrERGLOWLIBRARY 
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